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Abstract 

Multiplicity moments of charged particles in deep inelastic e~^p scattering have been 
measured with the ZEUS detector at HERA using an integrated luminosity of 38.4 
pb~^. The moments for > 1000 GeV^ were studied in the current region of 
the Breit frame. The evolution of the moments was investigated as a function of 
restricted regions in polar angle and, for the first time, both in the transverse mo- 
mentum and in absolute momentum of final-state particles. Analytic perturbative 
QCD predictions in conjunction with the hypothesis of Local Parton-Hadron Dual- 
ity (LPHD) reproduce the trends of the moments in polar-angle regions, although 
some discrepancies are observed. For the moments restricted either in transverse 
or absolute momentum, the analytic results combined with the LPHD hypothesis 
show considerable deviations from the measurements. The study indicates a large 
influence of the hadronisation stage on the multiplicity distributions in the restricted 
phase-space regions studied here, which is inconsistent with the expectations of the 
LPHD hypothesis. 
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1 Introduction 



The hadronic final state in deep inelastic scattering (DIS) events is the result of a hard 
partonic scattering initiated at large momentum transfers, Q 3> A, where A is a characteristic 
QCD scale of the order of a few hundred MeV. The subsequent parton cascade is followed by 
a soft fragmentation process. The latter occurs with small momentum transfers which may 
be considered to extend up to a value Qq, which is a QCD cut-off above which perturbative 
methods can be applied. There are two main approaches for describing the hadronic final 
state. The first comprises analytic perturbative QCD calculations based on the hypothesis 
of Local Parton-Hadron Duality (LPHD) this hypothesis states that parton-level QCD 
predictions are applicable for sufficiently inclusive hadronic observables, without additional 
assumptions about hadronisation processes. Therefore, the hadronic spectra are proportional 
to those of partons if the cut-off Qq is decreased towards a small value of A. The second 
approach is formulated in Monte Carlo (MC) programs, which generate a partonic final state 
according to the perturbative QCD picture and, below Qo — 1 GeV, hadronise the partons on 
the basis of non-perturbative models. The MC models are able to reproduce many detailed 
properties of the hadronic final state, but contain a large number of free parameters. 

Multiplicity distributions and correlations between final-state particles are an important 
testing ground for analytic perturbative QCD in conjunction with the LPHD hypothesis, 
as well as for MC models describing the hadronic final state 0. In this paper, a recent 
analysis of two-particle angular correlations is extended to a study of multiplicity dis- 
tributions measured in restricted phase-space regions of neutral current DIS events. The 
particle multiplicities are studied in terms of the normalised factorial moments^ 

Fg{n) = {n{n -l)...{n-q + g = 2, 3, . . . , 

for a specified phase-space region of size Q. The number, n, of particles is measured inside 
Q and angled brackets (...) denote the average over all events. The factorial moments, 
along with cumulants and bunching parameters ||^, are convenient tools to characterise 
the multiplicity distributions when Q becomes small. For uncorrelated particle production 
within Q, Poisson statistics holds and Fg = 1 for all q. Correlations between particles lead 
to a broadening of the multiplicity distribution and to dynamical fluctuations. In this case, 
the normalised factorial moments increase with decreasing Q. If the rise follows a power law, 
this effect is frequently called "intermittency" 0. 

The LPHD hypothesis, originally suggested for single-particle spectra, presently has experi- 
mental support in e'^e~ annihilations [0. However, the success of LPHD is less evident for the 
moments of inclusive single-particle densities at HERA energies P] and two-particle angular 
correlations 0. In contrast to the single-particle densities and global event-shape variables, 
the factorial moments are most directly affected by inclusive mani/-particle densities. This 
can be seen from the relation (n(n — 1) . . . (n — g + 1)) = J^p^'^\pi, . . .pqjdpi . . . dpq between 
the factorial moments in a region Vt and the inclusive g-particle density, p^'^^pi, . . .pg) = 
(l/cTtoi) d'^a/dpi . . . dpg. The latter represents the probability to detect q particles with 

^It may be noted that, for q = 1, Fq — 1 hy definitfon 
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momenta pi, . . .pg, irrespective of the presence of any other particles. Thus the factorial 
moments are an important tool to study such densities, for which the applicability of LPHD 
remains questionable even for higher e~^e~ energies at LEP [^[l^. Deviations of the mea- 
sured moments from theoretical predictions would imply the necessity for further refinements 
of perturbative QCD calculations and/or the failure of the LPHD hypothesis. In the lat- 
ter case, the magnitude of such deviations can shed light on details of the hadronisation 
processes. 



2 Analytic QCD results 

The dominant source of particle production inside a jet is gluon splitting in the QCD cascade, 
where the presence of a gluon enhances the probability for emission of another gluon nearby 
in momentum space. This leads to inter-parton correlations and non-Poissonian statistics for 
the multiplicity distributions in restricted phase-space intervals where partons are counted. 
Thus the factorial moments, Fg, are expected to deviate from unity. 

There exist a number of analytic QCD predictions for the moments of the parton multi- 
plicity distributions obtained in the Double Leading Log Approximation (DLLA), which are 
discussed below. 



2.1 Multiplicity-cut moments 

The normalised factorial moments of the multiplicity distributions of gluons which are re- 
stricted in either transverse momentum pt < p^"* or absolute momentum p = \ p |< p'^'^* are 
expected, respectively, to have the following behaviour |TT|: 



6 HE/Qo) ' 

(1) 

^ C{q)>l, 

where E is the initial energy of the outgoing quark that radiates the gluons and C{q) are 
constants depending on q. The maximum transverse momentum, p™*, and the pt values are 
defined with respect to the direction of this quark (see Fig. ^(a) for a schematic representation 
in the case of DIS). 

Equations (1) indicate that there are positive correlations between partons, because the 
factorial moments are larger than unity. For small p™* values, however, the correlations 
vanish due to the presence of an angular ordering of the partons in the jet. The existence 
of the angular ordering leads to the absence of branching processes with secondary gluon 
emission at small p^*^*. This ultimately leads to the suppression of the correlations in this 
phase-space region and to independent parton emission^ [Eg ^ 1 for p^^^ = Qq). On the 

^ Note that this effect is simiiar to muitiple photon brcmsstrahiung in QED. 
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other hand, the distribution of soft gluons with hmited absolute momenta, p < p^^^, remains 
non-Poissonian for any small value of p'^'^^. 

Note that the theoretical formulae (1) can be regarded as asymptotic calculations, i.e. they 
are valid at small p™* and p™*. Therefore, the calculations should be considered only as 
qualitative predictions when compared to the data using the LPHD hypothesis. 



2.2 Angular moments 

Gluon multiplicities in polar-angle rings around the outgoing-quark direction are expected 
to exhibit a distribution which is substantially wider than for Poisson statistics [^^-[T^. The 
moments rise with increasing values of the phase-space variable z as 

ln^ = ^a(l-/^.)(g-l), (2) 
where z and the constant a are defined as: 

-1, ©0 

z = a "^~Q^ 

(3) 

a = m — - — . 
A 

The angle 0o is the half opening-angle of a cone around the outgoing quark radiating the 
gluons, E is its energy, and A is an effective QCD scale (A 7^ A;^). A schematic represen- 
tation of the variables is shown in Fig. |l|(b), in which a window in the angular half- width 
6 is subtended by the ring centered at Gq. The factors Dq are the so-called Renyi dimen- 
sions |jl5[ known from non-linear physics. A decrease of the angular window corresponds 
to an increase of the variable z. The maximum possible phase-space region (0 = ©q), cor- 
responding to 2; = 0, is used to calculate the factorial moments, Fg{0), which are necessary 
for normalisation in Eq. (2). Such normalisation removes the theoretical uncertainty in the 
absolute values of Fq{z) and thus allows a comparison of the predicted ^-dependence of the 
moments with experimental results. 

For independent particle production in a restricted range of 26, Dq = 1 and Fq{z) = Fq{0) 
in Eq. (2). In contrast, the analytic QCD expectations for Dq obtained in the DLLA are as 
follows [ITM^ 



if the strong coupling constant, a^, is fixed at scale Q ~ EQq, then for large angular 
bins (i.e. small z), 

D, = 70^, (4) 



where 70 = y2CAas/vr is the anomalous QCD dimension, and Ca = 3 is the gluon 
colour factor. The as value is calculated at the one- loop level with n/ = 3 (number 
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of flavours). Such a behaviour of the moments corresponds to the expectations from 
intermittency, which can be interpreted in terms of the fractal structure of the gluon 
splittings p|. 



For running-a^ values, the Dq can be approximated by different expressions, depending on 
the treatment of the non-leading DLLA contributions to the moments: 



for Dokshitzer and Dremin [121: 



Dq ~ 7o 



q + l 



1 + 



4g2 ■ 



(5) 



for Brax et al. ITS 



q + l /I - 



(6) 



for Ochs and Wosiek 



q - w(q,z) 
z(g-l) ■ 



w 



{q,z) = gVl - z 1 



ln(l - z) 
2g2 



(7) 



an estimate for Dq which includes a correction from the Modified Leading Log Ap- 



proximation (MLLA) has also been obtained by Dokshitzer and Dremin [12|. In this 
case, Eq. (5) remains valid, but 70 is replaced by an effective jQ^iq) which depends on 

q- 



All the theoretical calculations quoted above are performed for partons. Therefore, for 
comparisons of the predictions with data, LPHD is assumed. 



3 Experimental set-up 

The data were taken in 1996 and 1997 using the ZEUS detector at HERA. During this 
period, the energy of the positron beam, E^, was 27.5 GeV and that of the proton beam was 
820 GeV. The integrated luminosity used for the present study was 38.4 ± 0.6 pb"'^. 

ZEUS is a multi-purpose detector described in detail elsewhere [|T6|. Of particular impor- 
tance for the present study are the central tracking detector (CTD), positioned in a 1.43 T 
solenoidal magnetic field, and the calorimeter. 

The CTD is a cylindrical drift chamber with nine superlayers covering the polar angle^ 
region 15" < 9 < 164° and the radial range 18.2-79.4 cm. Each superlayer consists of eight 

^ The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, and the X axis pointing left towards the centre of HERA. The coordinate origin is 
at the nominal interaction point. The polar angle 9 is defined with respect to the positive Z-direction. The 
pseudorapidity, 77, is given by — ln(tan |). 
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anodes. The transverse-momentum resolution for charged tracks traversing all CTD layers is 
= 0.0058p^ © 0.0065 © 0.0014/p_|_, with being the track-transverse momentum 
(in GeV). The single hit efficiency of the CTD is greater than 95%. 

The CTD is surrounded by the uranium-scintillator calorimeter (CAL), which is longitudi- 
nally segmented into electromagnetic and hadronic sections. The relative energy resolution 
of the calorimeter under test beam conditions is 0.18/\/^ for electrons and 0.35/a/E' for 
hadrons (with E in GeV). The CAL cells provide time measurements with a time resolution 
below 1 ns for energy deposits larger than 4.5 GeV. 



4 Data selection 



4.1 Event kinematics 

The basic event kinematics of DIS processes are determined by the negative squared four- 
momentum transfer at the positron vertex, = —q^ = —{k — k ^ {k and k denote the 
four-momenta of the initial- and final-state positrons, respectively) and the Bjorken scaling 
variable x = Q^/(2P • q), where P is the four-momentum of the proton. The fraction of 
the energy transfered to the proton in its rest frame, y, is related to these two variables by 
y ^ Q^/ sx, where ^/s is the positron-proton centre-of-mass energy. 



The moments were studied in the Breit frame |T^, which provides a natural system to 



separate the radiation from the outgoing struck quark from the proton remnants. In this 
frame, the exchanged virtual boson with virtuality is completely space-like and has a 
momentum q = (q'o, ^Xb? ^Ib; '?^b) — (0)0)0)~Q)- the quark-parton model, the incident 
quark has pzj^ = Qjl and the outgoing struck quark carries "pz^ = — Q/2. All particles with 
negative pzj^ form the current region. These particles are produced by the fragmentation 
of the struck quark, so that the current region is analogous to a single hemisphere in e~^e~ 



annihilations. The variables pt = \JPx^ ~^ Pyb ^"^^ ®o were determined with respect to the 
negative Ze-direction, which, for high events, is a good estimate of the outgoing struck- 
quark direction. The moments in the current region of the Breit frame were measured in the 
kinematic region > 1000 GeV^, which is sufficiently high to avoid major contributions 
from the boson-gluon fusion (BGF) processes (see Section 0). 



4.2 Event and track selection 

The kinematic variables x and were determined by a combination of methods: (i) the 
measurement of the energy and angle of the scattered positron (denoted by the subscript 



e); (ii) the double angle (DA) method [0 and (iii) the Jacquet-Blondel (JB) method ||T! 
depending on which is most appropriate for a particular cut. In the DA method, the variables 
X, and y were reconstructed using the angles of the scattered positron and the hadronic 
energy flow. In the JB method, used only in the event selection, the kinematic variables 
were determined entirely from the hadronic system. 



5 



The boost vector from the laboratory to the Breit frame was determined using the DA 
method, which is less sensitive to systematic uncertainties in the energy measurement than 
the other methods. In the reconstruction of the Breit frame, all charged particles were 
assumed to have the pion mass. 

The triggering and online event selections were identical to those used in a previous publi- 



cation 1 20]. To preselect neutral current DIS events, the following cuts were applied: 



• E^>10 GeV, being the energy of the scattered positron; 

• QIa > 1000 GeV2; 

• 35 < (5 = "^Ei^l — cos 9i) < 60 GeV, where Ei is the energy of the i*^ calorimeter 
cell and 6i is its polar angle with respect to the beam axis; 

• Ve < 0.95; 

• ym > 0.04; 

• Z coordinate of the reconstructed vertex position, as determined from the tracks fitted 
to the vertex, in the range —40 < Z^^rtex < 50 cm; 

• a timing cut requiring that the event time as measured by the GAL be consistent with 
an e~^p interaction. 



The resulting data set contained 7369 events. 



The present study is based on the tracks reconstructed with the GTD and fitted to the 
event vertex. The scattered positron was removed from the track sample. In addition, the 
following cuts to select tracks were imposed: 

• p'f''^ > 0.15 GeV, where p'f''' is the transverse momentum of the track; 

• \r]\< 1-75. 

These cuts restrict the data to a well understood, high-acceptance region of the GTD. 



5 Event simulation 



Neutral current DIS events were generated with ARIADNE 4.10 with the high 
modification developed during the HERA Monte Garlo workshop using tuned parameters 
23i . This MG program is based on the colour-dipole model, in which gluons are emitted from 



the colour field between quark-antiquark pairs, supplemented with the BGF process. The 



hadronisation was simulated using the Lund string model as implemented in JETSET p4 
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Charged hadrons with a mean hfetime longer than 10"*^ seconds were considered as stable. 
The CTEQ4D |^ parameterisation of the proton parton distribution functions was used. 
The events were generated without any cuts on kinematic variables and track quantities. 
The MC events obtained by this procedure defined the generator-level sample. 

To obtain a detector-level sample, events were generated with the DJANGO 6.24 pro- 



gram based on HERACLES 4.5.2 ^7j, in order to incorporate first-order electroweak cor- 
rections. The parton cascade and hadronisation were then simulated with ARIADNE. The 



events were then processed through a simulation of the detector using GEANT 3.13 |2^ to 
take into account particle interactions with the detector material, particle decays, event and 
track migrations, resolution, acceptance of the detector and event selection. The detector- 
level MC events were processed with the same reconstruction program as was used for the 
data. 



In addition to ARIADNE, the LEPTO 6.5 [||] and HERWIG 5.9 |30| generators with tuned 



parameters were also used to compare with the data. The parton cascade in LEPTO is 
based on a matrix-element calculation matched to parton showers according to the DGLAP 
equations. As in ARIADNE, the hadronisation in LEPTO is simulated with JETSET. 
HERWIG also has a parton shower based on DGLAP evolution, but parton emissions are 
ordered in angle, in contrast to LEPTO, which orders the parton emissions in invariant mass 
with an additional angular constraint to ensure coherence. The hadronisation in HERWIG 



is described by the cluster fragmentation model |^ 



The Bose- Einstein (BE) interference between identical particles can produce additional cor- 
relations which distort the behaviour of the moments measured in small phase-space regions. 
By default, the BE interference is turned off in JETSET; it is absent in HERWIG. The effect 
of the BE correlations on the present results is considered in the next section. 



6 Correction procedure 

The bin- by-bin correction factors, Cg = Af^^/A^^^, for a given kinematic region in were 
evaluated using the ARIADNE 4.10 MC event samples, separately for each observable, 
Aq = and Fg{z)/Fq{0), where ^f" (^^*^*) was calculated at the generator 

(detector) level. The factors Cg correct the data for the detector effects and event-selection 
cuts as well as for the effects that are not included in the analytic calculations, namely (i) 
Dalitz decays of the vr"; (ii) leptons in the final state and (iii) initial-state photon radiation. 

The correction factors for the p™* and p*^"* moments are approximately fiat and vary between 
about 1.0 and 1.3, except for the smallest momentum cuts and the highest-order p™* and 
pcut ]2ioments, where the correction factors are as large as 1.7. The correction factors for the 
Fg{z)/Fg{0) moments vary smoothly from 1.0 to 1.5. The contribution of Dalitz pairs to the 
correction factors is negligible, except for the smallest phase-space regions in p™*, p'^'^* and 
0, where it does not exceed 20%. 

The analytic QCD calculations do not include BE interference. To take into account this 
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phenomenon, the correction factor, BEg, was calculated for each data point 



■on 



= BE, gen ' (8) 

where Af^ {A^^'^^^) are the moments at the generator level without (with) the BE effect. 
For this correction, the BE effect according to a Gaussian parameterisation with JETSET 
default parameters was used. As a cross check, the BE correction was determined from 



the ZEUS data and also using the parameters obtained by the HI Collaboration The 
differences between the correction factors BE, for these three sets of parameters were small. 
The values of the BE correction factors are discussed in Section |8|. 



7 Systematic uncertainties 

The overall systematic uncertainty for each moment was determined from the following 
uncertainties added in quadrature: 

• the choice of the MC models. HERWIG was used in place of ARIADNE to determine 
the correction factors. This uncertainty was typically 45% of the total systematic 
uncertainty; 

• event reconstruction and selection. Systematic checks were performed by varying the 
cuts on S, Ue, yjB and the vertex position requirement (40 < 6 < 55 GeV, ye < 0.85, 
yjB > 0.05 and —35 < ^vertex < 45 cm). The overall uncertainties associated with the 
event selection were typically 25% of the total systematic uncertainty; 

• track reconstruction and selection. Tracks were required to have transverse momenta 
larger than 0.2 GeV and \ f] \< 1.5. Tracks that reach at least the third CTD superlayer 
were used. These uncertainties were typically 30% of the total systematic uncertainty. 

As an example, the total systematic uncertainty of F2{p™^) for the smallest p^*^* value used in 
this analysis is approximately 6%. For the large pj:"* regions, a typical value of the systematic 
uncertainty is 1%. 

To ensure that the two particles in a pair were resolved, the angle between two tracks in 
the current region of the Breit frame should not be smaller than two degrees 0]. This limit 
determines the minimum size of phase-space regions used in this analysis. 

The kinematic variables pt, and Gq have to be measured with respect to the outgoing- 
quark direction, which coincides with the ZB-direction of the Breit frame only for the zeroth- 
order 7*g scattering. The first-order QCD effects, BGF {'y*g —>■ qq) and QCD Compton 
~^ Q9)i can lead to two hard partons that are not coUinear to the ZB-direction. For 
these events, the determination of the kinematic variables with respect to the Z-q axis of 
the Breit frame is inappropriate, especially when these two hard partons escape the current 
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region. Using the LEPTO MC model, it has been determined that, at > 1000 GeV^, 
the number of BGF and QCD Compton events without hard partons in the current region, 
relative to the total number of events, is < 15%. As a systematic check, the factorial moments 
have been calculated in kinematic variables with respect to the thrust axis determined from 
the charged particles in the current region of the Breit frame. The results agree with those 
for which the negative Ze-direction of the Breit frame is used. 



8 Results 

The charged-particle multiplicities measured in the entire current region without additional 
momentum cuts on the charged particles were compared to the previous ZEUS measurements 
and to e~^e~ data |Q. For the present measurement, the corrected average value of is 



(Q^) — 2070 GeV^, which corresponds to an energy v/(Q2) ~ 46 GeV. The corrected average 
charged-particle multiplicity in the current region, 5.76 ± 0.29, where the error represents 
statistical and systematic uncertainties added in quadrature, agrees with the previous ZEUS 
measurement, 6.01 ± 0.46, for a similar value of (Q^). This average multiplicity is lower 
than that found in a single hemisphere of e~^e~ at ^/s = 44 GeV, which is 7.58 ± 0.37 
33| . The second-order factorial moment F2 in a single hemisphere of e~^e~ annihilations isQ 



F2 = 1.036 ±0.026, which is consistent with the value 1.051 ±0.018 obtained in this study of 
the entire current region. This difference in the average multiplicity can be attributed to a 
migration of particles between the current and target regions due to BGF and QCD Compton 
processes, an effect which has been theoretically quantified and directly measured by 



ZEUS 0]. Such a migration leads to current-target anti-correlations which do not vanish 
even at the high values studied here. This is unlike the e~^e~ annihilation processes, where 
weak and positive forward-backward correlations have been observed at a similar energy BH] . 



8.1 Multiplicity-cut moments 

The multiplicity moments of order q = 2, . . . , 5 as a function of p1^^ are shown in Fig. |^. The 
pcut ]2ioments were calculated with respect to the negative Ze-direction of the Breit frame. 
As pj"* decreases below 1 GeV, the moments rise. This disagrees with the DLLA theoretical 
calculations, which predict that the moments approach unity (see Eq. (1)). On the other 
hand, Monte Carlo models show the same trend as the data, although some differences in 
slopes are apparent. LEPTO exhibits smaller values of the moments than for the data at 
small pI'^^ values, while HERWIG overestimates the moments. ARIADNE gives the best 
overall description of the correlations. 

Figure also shows the factors BEg from Eq. (8) used to correct the data for the BE 
interference. The largest correction was found in the case of moments of high orders (g = 4, 5) 
for small p™^ cuts (< 0.5 GeV). 

The second-order moment in a single hemisphere of e+e" annihilation can be estimated from the 
dispersion D — y/ {n?) — (n)^ and the average multiplicity {n) using the relation F2 = 1 + ("ra)"^ — 



9 



To investigate contributions of the hadronisation processes to the moments, further sam- 
ples of the ARIADNE MC were generated: (1) considering neutral hadrons in addition to 
charged particles; (2) using only hadrons coming from the Lund string fragmentation with- 
out resonance decays. For these two cases, the moments have an even steeper rise at low 
j9^"*. Replacing the Lund string- fragmentation scheme by independent fragmentation and 
removing resonance decays led to an increase of the moments for p^"* ~ 0.5 GeV, but below 
this cut the moments approach constant values which are larger than unity. An MC sample 
generated without the colour-dipole parton emission and resonance decays, but retaining the 
Lund string fragmentation, leads to an increase of the moments with decreasing p^^^. 

Figure ^ shows the factorial moments for restricted p values. As was the case for the p™* 
moments, there is no indication that the values approach unity at small p'^"* values. All MC 
models at the hadron level agree with the data. The BE interference has a large influence 
on this measurement, giving a 20% correction. 

The DLLA QCD calculations neglect energy-momentum conservation in gluon splittings, 
considering idealised jets at asymptotically large energy. This effect can be taken into ac- 
count in the Modified Leading Log Approximation; however, such calculations have not yet 
been performed. Therefore, it is important to compare the experimental results with the 
parton shower of MC models that explicitly include energy-momentum conservation and the 
complete parton-splitting functions. In addition, such comparisons allow the investigation 
of the contributions from hadronisation. 

The moments restricted in pt and p variables were determined from the ARIADNE parton- 



level, the physics implementation of which strongly resembles the analytic calculations ||TT 
The LPHD requires Qo to approach A. The ARIADNE defauh value of A is 0.22 GeV, while 
the transverse momentum cut-off Qo is 0.6 GeV. Therefore, to satisfy the LPHD requirement, 
the parton predictions were obtained after switching off the string fragmentation and lowering 
the Qo value in ARIADNE to Qo = 0.27 GeV. This value was chosen in order to insure that 
the average multiplicity of partons in the current region of the Breit frame equalled that of 
hadrons, which is important for normalisation purposes. This method is somewhat different 
to that used in a similar study of the ARIADNE parton cascade for e~^e~ annihilation 



events [0]. In contrast to that approach, the present modification of the model retains the 
original default value of A, to preserve the kinematics of first-order QCD processes in the 
Breit frame. 

The parton-level ARIADNE for p^*^* and p'^"* moments is shown as thin solid lines in Figs. ^ 
and Agreement between the parton-level predictions and the data for the lowest pj"* 
values can be excluded at greater than 99% confidence level for q = 2,3 and at greater 
than 95% confidence level for q = 4,5. The data for all p™* moments disagrees with the 
parton-level predictions by many standard deviations. The MC predictions for partons have 
the trends expected from the analytic DLLA results, showing a decrease of the moments 
for decreasing p^*^* and p'^^^. This is contrary both to the measurements and to the hadron- 
level MC predictions, which show a rise of the moments. Analogous differences between the 



hadron and parton levels of ARIADNE have also been observed in e^e annihilation ITJ] 



The parton predictions in the full current region (large p^"* and p'^^^^ show larger values of 
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the moments than the data, despite the tuning of the parton to the hadron multiphcities. In 
the case of e"'"e~, the MC-tuned parton level describes the hadron- level predictions for large 

It is likely that the discrepancy between partons and hadrons in DIS 



momentum cuts 11 



is due to additional non-perturbative effects related to the proton remnant. 



8.2 Angular moments 

Figure ^ shows the angular moments Fq{z) / Fq{Qi) . This represents the first measurement 
of these quantities in DIS. The measurements are compared to DLLA QCD calculations in 
Eq. (2), after substitution of different forms of the Renyi dimension (see Eqs. (4)- (7) ) and 
using the effective QCD parameter A = 0.15 GeV as in LEP studies For DIS in the 

Breit frame, the energy E = Q/2 was chosen to define the ag value, corresponding to the 
outgoing quark momentum in the quark-parton model. The factors used to correct the data 
for the BE interference are close to unity for all angular moments (not shown). A significant 
disagreement with the data is observed, which becomes smaller for higher orders of the 
moments. Although the error bars are clearly correlated, even for the high-order moments, 
agreement between the predictions and the data can be excluded for z =0.05-0.2 at greater 
than 95% confidence level. For running-as calculations, the best description is observed for 
high-order moments in the case of the DLLA result with a correction from the MLLA. This 
result is similar to that obtained for e^e~ annihilation It should be noted that the 

present analysis, however, uses a larger angle Gq = 45°, in order to increase the statistics, 
and that the average energy {{E) ~ 23 GeV) of the outgoing quark giving rise to the parton 
shower is smaller than at LEP by a factor of two. 

Figure |^ shows the comparisons of the angular moments with Monte Carlo models at the 
hadron level. All models reproduce the trends of the data but somewhat overestimate the 
strength of the correlations; however, LEPTO fails to give a reasonable description of the 
angular moments. 

The parton-level ARIADNE predictions for the angular moments are shown in Fig. |^. The 
ratio Fq{z) / Fq{{]) predicted by ARIADNE at the parton level (thin lines) overestimates the 
data, but is closer than the analytic DLLA calculations for the second-order moment. This 
comparison with the MC parton level and similar studies of e^e~ annihilations at LEP [p|, p!0|] 
imply that MLLA analytic calculations are necessary to give a better description of the 
measured angular moments. 



9 Conclusions 

Multiplicity moments have been studied in deep inelastic scattering at > 1000 GeV^ 
with the ZEUS detector at HERA. The moments measured in the current region of the Breit 
frame show a strong rise as the size of the restricted intervals in pt, p and G decreases. Monte 
Carlo models, which include hadronisation effects, reproduce the correlation pattern of the 
hadronic final state, although there are some discrepancies with the data. 
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The analytic QCD calculations for the angular moments agree qualitatively with the data, 
but there are discrepancies at the quantitative level. An inadequate treatment of energy- 
momentum conservation in gluon splittings, which is important at low energies, is likely 
to be responsible for the discrepancies between the analytic QCD expectations and the 
data. The parton level of ARIADNE, which does not suffer from the limitations of the 
DLLA calculations, indicates that the discrepancies can be reduced. MLLA calculations are 
required to provide an accurate description of the angular moments. Similar studies of the 
angular moments measured at LEP also show that the analytic perturbative calculations 
need to be improved. 

In contrast to the angular fluctuations, the moments calculated in restricted pt and p regions 
show significant discrepancies, both for the DLLA calculations and for the parton-level of 
ARIADNE. This is most clearly seen for thep^"* moments; while the analytic calculations and 
ARIADNE predictions for partons decrease and approach unity for small pt regions, the data 
show a steep rise. This is the first indication that perturbative QCD, in conjunction with the 
LPHD hypothesis, fails on a qualitative level to describe the hadronic multiphcities. Monte 
Carlo studies of the hadronisation effects indicate that this is due to violation of the LPHD 
hypothesis for many-particle densities measured in the regions restricted in the transverse 
momenta with respect to the direction of the outgoing quark. 
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a) 




Figure 1: A schematic representation of the measurements of the factorial moments in DIS: 
(a) restricted in the transverse momentum, pt, defined with respect to the outgoing-quark 
direction and absolute momentum, | p |; (b) in polar-angle rings of width 20. 
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Figure 2: Factorial moments for charged particles in the current region of the Breit frame 
as a function of p™^, compared to Monte Carlo models. The thin lines show the parton-level 
prediction of ARIADNE with Qq = 0.27 GcV and the thick solid lines are the hadron-level 
predictions of the same model. The dashed and dotted lines show the hadron-level predictions 
of HERWIG and LEPTO, respectively. The inner error bars are statistical uncertainties; the 
outer are statistical and systematic uncertainties added in quadrature. The correction factors 
for the BE effect are shown below each plot. 
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Figure 3: Factorial moments for charged particles in the current region of the Breit frame 
as a function of p™^, compared to Monte Carlo models. The thin solid lines show the parton- 
level prediction of ARIADNE and the thick solid lines are the hadron-level predictions of the 
same model. The dashed and dotted lines show the hadron-level predictions of HERWIG 
and LEPTO, respectively. The inner error bars are statistical uncertainties; the outer are 
statistical and systematic uncertainties added in quadrature. The correction factors for the 
BE effect are shown below each plot. 
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ZEUS 




Figure 4: Fq(z)/Fq(0) (g = 2, . . . , 5) as a function of the scaling variable z for 6o = 45° and 
A = 0.15 GeV, compared to the QCD prediction of Eq. (2) with different parameterisations 
of Dq. as = const of Eq. (4) (bold sohd lines); (1) DLLA of Eq. (5) (dashed lines); (2) 
DLLA of Eq. (6) (dotted lines); (3) DLLA of Eq. (7) (dash-dotted lines); and (4) DLLA 
with a correction from the MLLA |jl2|] (thin solid lines). The values of B are given at the 
top of the plots. The inner error bars are statistical uncertainties; the outer are statistical 
and systematic uncertainties added in quadrature. 
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Figure 5: The data of Fig. ^ compared to different Monte Carlo models at the hadron level: 
ARIADNE (bold solid lines), HERWIG (bold dashed lines) and LEPTO (bold dotted lines). 
The thin solid lines show the parton-level prediction of ARIADNE with Qq = 0.27 GeV. 
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